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We describe a non-perturbative determination of ca using correlators containing the axial-vector and pseu- 
doscalar currents at zero and non-zero momentum. We apply the method of Bhattacharya et al to extract ca 
from the requirement that the ratio of appropriate correlators for the PCAC relation becomes independent of 
time in the excited state region. We find that the result depends strongly on the order of the derivatives used 
in the PCAC relation. We also find that, using the lowest order derivatives, we cannot get a consistent value of 
ca between zero and non-zero momentum cases. The ca values that we obtain as we improve the derivatives are 
consistent and decrease in magnitude heading towards the perturbative result. 



1. Introduction 

The use of Symanzik-improved lattice actions 
and matrix elements is widespread. However, 
with each improvement term added the corre- 
sponding coefficient must be determined to enable 
discretisation effects to be reduced to the desired 
level. Considering the light hadron spectrum and 
matrix elements, the relevant 0(a) improvement 
coefficients are, for the most part, only known to 
1-loop in perturbation theory. A nonperturbative 
determination of these coefficients is desirable and 
may be simpler than performing higher loop per- 
turbative calculations. 

The ALPHA collaboration, using Schrodinger 
functional techniques, have calculated several 
0(a) improvement coefficients nonperturbatively. 
In most cases agreement is found with 1-loop 
(tadpole-improved) perturbation theory, or the 
discrepancy is consistent with estimates of the 
size of omitted higher orders in a (albeit in some 
cases assuming a slow convergence of the pertur- 
bative series). However, for the 0(a) improve- 
ment of the axial-vector current 
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where = "07^ 75^ and P = ipjsip, ca is found 
to be many times larger than the 1-loop pertur- 
bative value even at reasonably high 0, as shown 
in figure |l|, taken from reference |l| . 

ca appears in the expression for the pion decay 




Figure 1. Ca extracted by the ALPHA collabo- 
ration. The circles indicate the results of refer- 
ence 0, the dotted line the results of 1-loop per- 
turbation theory and the crosses the results from 
1-loop tadpole-improved perturbation theory. 



constant, f n , 

U = (1 + b A )Z A (fi 0) + c A afM + 0(a 2 ) (2) 

and the bare quark mass (from the PCAC rela- 
tion) 
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Choosing the perturbative or nonperturbative 
value for c A leads to significantly different val- 
ues for these quantities at finite (3. Bhattacharya 
et. al. proposed an alternative method for deter- 
mining ca, and other 0(a) coefficients in refer- 
ence Here, we apply their method, but with 
higher statistics than reported in reference M . In 
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addition, we extract ca at finite, as well as zero, 
momentum. 

2. C A from the PCAC relation 

The PCAC relation in euclidean space 

< d^A^P >= 2m < P(x)P > (4) 

should hold for all x on the lattice up to discreti- 
sation terms. For simplicity we assign 



r(f) = 

S(t) : 

Thus, 
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< dlP{t)P > 



< P(t)P > 



(5) 
(G) 



2m = r(t) + 0(a) 

2m = r{t) +c A s(t) + 0(a 2 ). 



(7) 
(8) 

Equation || only holds if the 0(a)-improvement 
term (the clover term) is included in the light 
quark action and the coefficient, c$Wi is deter- 
mined nonperturbatively. The method of Bhat- 
tacharya et. al. is to determine ca by minimizing 
the dependence of r(t) + CAs(t) on t and, hence, 
reduce the discretisation errors in this quantity. 
Using this method, we investigate the effect of im- 
proving on the determination of ca] discreti- 
sation errors in the lattice derivatives may be the 
dominant source of lattice spacing errors in the 
bare quark mass, and hence, the value obtained 
for ca may be artificially high. 

Normally, the symmetric lattice derivatives 
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are used, which contain 0(a 2 ) errors. We will 
consider the effect of using 



A(+-) = a<+-> - ^A+A(+- 
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which are correct up to C*(a 4 ). One can continue 
to correct to 0(a 6 ), we denote these derivatives 

A' + and A'„ . As an additional constraint on 
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Table 1 

Simulation details. 



,8 = 6.0 

o LL: A 0(a z ) 
FL: A 0(a z ) 
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Figure 2. The ratio of correlators r(t). 



ca, we investigate whether Lorentz invariance is 
restored to sufficient accuracy i.e. that consis- 
tent values of m are obtained at zero and finite 
momentum. 

3. Estimating Ca using excited states 

The configurations and light quark propagators 
were provide by the UKQCD collaboration. The 
simulation parameters are given in table |l|. For 
more details see reference |J . Note that the Csw 
values correspond to those determined by the AL- 
PHA collaboration. 

Figure || shows the ratio of correlators, r(t), 
for the data set at j3 = 6.0 for correlators local 
at the source and sink (LL) and those fuzzed at 
the source and local at the sink (FL). The 0(a 2 ) 
lattice derivatives, equation [| were used. To- 
wards the center of the lattice r(t) tends to a 
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Figure 3. The 0(a)-improved bare quark mass as 
determined using from LL correlators. 



constant, 2m ~ 0.1096. However, in the LL case, 
close to the origin there are significant deviations 
from this value due to discretisation effects. The 
much wider plateau in r(t) for the FL correla- 
tors (where the fuzzing has been chosen to more 
or less eliminate excited state contributions) show 
that the discretisation effects in r(t) LL are asso- 
ciated with excited states. We can determine ca 
by adding the term CAs(t) to reduce the discreti- 
sation effects at small t. To do this we perform 
a correlated fit to r(t) using r(t) — 2m — CAs(t), 
to extract 2m and ca- The values of ca obtained 
are detailed in the next section. 

To see how well the discretisation effects are re- 
duced we plug these values back into r{t) +CAs(t) 
for LL correlators and compare with the unim- 
proved case. At (3 = 6.0 we find that the plateau 
in 2m can be extended from i = lltot = 6 
when ca = —0.10 and 0(a 2 ) lattice derivatives 
are used. However, as shown in figure ||, if one 
improves <9 M , the discretisation effects can be re- 
moved by the same amount but lower values of ca 
are required. Clearly the value of ca and the form 
of the derivatives significantly affects the value of 
m extracted. Furthermore, we find that consis- 
tent values of m at zero and finite momentum can 
only be obtained if 0(a A ) or 0(a 6 ) derivatives are 
used. Figures f| and || show the difference in the 
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Figure 4. The difference in the bare quark mass 
determined at zero and finite momentum, where 
standard lattice derivatives are used. 




Figure 5. The difference in the bare quark mass 
determined at zero and finite momentum, where 
improved lattice derivatives are used. 
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bare quark mass at zero and finite momentum for 
standard and improved derivatives. 

4. Results for Ca 

Results are shown in figures |and§. The main 
points are: 

• A reliable estimate for ca is indicated by 
a stable value with t m i n . At (3 = 6.0, the 
statistical errors grow rapidly for t m i n > 8, 
although the results are reasonably stable 
for t m i n = 6 — 8. The situation is similar at 
(3 = 6.2. 

• Clearly ca decreases when the temporal 
derivatives are improved. The decrease is 
less going from derivatives correct to 0(a 4 ) 
to those correct to 0(a 6 ), as expected. 

• ca obtained using improved derivatives is 
below the value obtained by the ALPHA 
collaboration. For the (9(a 6 )-correct deriva- 
tives ca is close to the 1-loop tadpole- 
improved value (although the perturbation 
theory has not been performed with im- 
proved derivatives). In principle, our re- 
sults can differ with those of reference JlJ] 
by 0{A QCD a) ~ .300/2.1 = .14 at f3 = 6.0 
and .10 at = 6.2. 

• The results can be compared to those ob- 
tained by Bhattacharya et. al. 0. They 
obtain c A = -0.02(2) (2) at = 6.0 and 
-0.033(4) (3) at = 6.2 for k c with a sig- 
nificantly smaller number of configurations. 
Our results are at finite k around n strange- 
A preliminary study indicates that the light 
quark mass dependence of ca is small. 

5. Conclusions 

We have applied the method of reference 
to extract the 0(a) improvement coefficient, ca- 
We find that the accuracy and reliability of this 
method is limited by the small range of timeslices 
over which a stable value of ca is found. Never- 
theless, we have shown that the value of ca ob- 
tained depends significantly on the discretisation 
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Figure 6. The results for ca at (3 = 6.0 obtain 
using correlated fits with Q > .10 (the dashed 
points have 0.10 > Q > .01). The data has 
been averaged over positive and negative (i.e. 
t > 25) timeslices. The errors are bootstrapped 
over 100 bootstrap samples. The result of the AL- 
PHA collaboration is shown as a horizontal line, 
where the dashed lines indicate the numerical er- 
ror. The tadpole-improved 1-loop perturbative 
result is also indicated, the error is taken to be 
1q5(tt/o). 



chosen for the lattice derivatives appearing in the 
PCAC relation, ca reduces as the lattice deriva- 
tives are improved, and lies between the value 
obtained by the ALPHA collaboration and the 1- 
loop tadpole-improved value. Furthermore, con- 
sistent values for the bare quark mass at zero and 
finite momentum can only be obtained for deriva- 
tives correct to 0(a 4 ) and above. 

These results suggest a similar study should 
be performed within the Schrodinger functional 
approach. In the future, we plan to investigate 
the quark mass dependence of ca and check the 
scaling properties of the renormaliscd light quark 
mass and f n obtained using our values of ca- 
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Figure 7. The results for c A at ft = 6.2. The 
figure is labelled in the same way as figure ^ 
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